Relaxation Scenarios in a Mixture of Large and Small Spheres: 
Dependence on the Size Disparity 
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We present a computational investigation on the slow dynamics of a mixture of large and small 
soft spheres. By varying the size disparity at a moderate fixed composition different relaxation 
scenarios are observed for the small particles. For small disparity density-density correlators exhibit 
moderate stretching. Only small quantitative differences are observed between dynamic features 
for large and small particles. On the contrary, large disparity induces a clear time scale separation 
between the large and the small particles. Density-density correlators for the small particles be- 
come extremely stretched, and display logarithmic relaxation by properly tuning the temperature 
or the wavevector. Self-correlators decay much faster than density-density correlators. For very 
large size disparity, a complete separation between self- and collective dynamics is observed for the 
small particles. Self-correlators decay to zero at temperatures where density-density correlations are 
frozen. The dynamic picture obtained by varying the size disparity resembles features associated 
to Mode Couphng transition lines of the types B and A at, respectively, small and very large size 
disparity. Both lines might merge, at some intermediate disparity, at a higher-order point, to which 
logarithmic relaxation would be associated. This picture resembles predictions of a recent Mode 
Coupling Theory for fluids confined in matrixes with interconnected voids [V. Krakoviack, Phys. 
Rev. Lett. 94, 065703 (2005)]. 

PACS numbers: 82.70.Dd, 64.70.Pf, 83.10.Rs 



I. INTRODUCTION 

Relaxation dynamics in glass-forming systems can 
be strongly sped up or slowed down by the addition of 
a component of very different mobility. An example 
is provided by mixtures of large and small particles 
of very different size. Several investigations have evi- 
denced the possibility of melting a glass of large parti- 
cles by adding a proper amount of small particlesiiSiS. 
Less attention has been paid to the dynamic fea- 
tures exhibited by the small particles. Very recently, 
we have presented a molecular dynamics investiga- 
tion on the relaxation dynamics of a mixture of large 
and small soft spheres^. We define the "size dispar- 
ity" , 5, as the ratio of the diameters of the large and 
the small particles. For a sufficiently large disparity 
[5 = 2.5 in Ref»^), at moderate and low concentrations 
of the small particles, the latters exhibit unusual re- 
laxation features. Differently from the usual two-step 
increase and decay observed for, respectively, mean 
squared displacements and dynamic correlators^iii^, 
these quantities do not exhibit a defined plateau at 
intermediate times between the microscopic and dif- 
fusive regimes^. This result suggests a softer char- 
acter for the collective caging mechanism — i.e., the 
temporary trapping of each particle by its neighbors. 
Dynamic correlators show a highly stretched decay, 
and for selected values of the control parameters the 
decay is logarithmic in time. By varying wavevec- 
tors or control parameters the decay shows a striking 
concave-to-convex crossover—. 

These anomalous relaxation features resem- 
ble predictions of the Mode Coupling Theory 
(MCT)ii2.i^!^iii for state points close to higher-order 
MCT transitions, initially predicted by schematic 
models^^, and later derived for one-component 
systems of particles interacting through a re- 



pulsive potential complemented by a short-ranged 
attractioni^ii^ii^ii^iii. These models are used as effec- 
tive pictures for the large particles. The short-ranged 
attraction arises as an effective interaction between 
the large particles, originating from the depletion 
mechanism induced by the presence of the small 
particlesi^i^. In a certain range of density, temper- 
ature and mixture composition, competition occurs 
between two different mechanisms for dynamic arrest 
of the large particles: steric repulsion characteristic of 
colloidal systems, and formation of reversible bonds, 
induced by the effective short-ranged attraction. The 
higher-order MCT scenario arises as a consequence 
of these two competing mechanisms of very different 
localization lengths^^'^^. When heating up or cooling 
down the system, dynamic arrest is exclusively 
driven by, respectively, steric repulsion and reversible 
bond formation, and relaxation features of standard 
liquid-glass transitions are recovered22i2i. Many of 
the predictions associated with the higher-order MCT 
scenario for short-ranged attractive colloids has been 
succesfuUy tested by simulation o^^i^^i^^i^^'^'^i^^ and 
experiments2^^'21^'29^. 

If the size disparity is sufficiently large to induce a 
large time scale separation between large and small 
particles, mean squared displacements and density- 
density correlators for the small particles in the mix- 
ture of Ref.^ display striking similarities with quali- 
tative features associated to higher-order MCT tran- 
sitions. Similar results have been reported for the 
fast component in simulations of polymer blends with 
components of very different mobilities^i^. These 
analogies with relaxation features in short-ranged at- 
tractive colloids suggest that the higher-order MCT 
scenario might be a general feature of systems show- 
ing a competition between different mechanisms for 
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dynamic arrest. For the case of the small particles 
in the mentioned polymeric and non-polymeric mix- 
tures, we have suggested a competiton between bulk- 
like caging and confinement^!^. Bulk-like caging is 
induced by neighboring small particles and confine- 
ment is induced by the slow matrix formed by the 
large particles. 

In this article we present new support for this in- 
terpretation by evidencing dynamic analogies with re- 
cent MCT theoretical calculations by Krakoviack-^'-'^ 
in a mixture of mobile and fixed particles, which cx- 
plicitely shows a higher-order MCT transition. We 
have carried out new simulations, at a fixed mixture 
composition, for disparities 1 6 < 8, extending re- 
sults reported in Ref.^ for 6 = 2.5. For 6^1 standard 
relaxation features are recovered. For the largest in- 
vestigated values of (5 a new relaxation scenario arises 
for the small particles, showing features characteristic 
of nearby MCT transitions of the so called type-A. 
Such transitions are defined by a zero value of the 
long-time limit (critical non-ergodicity parameter) of 
density-density correlators, different from the finite 
value (observed at smaller S) defining the usual type-B 
transitions. This feature agrees with early MCT pre- 
dictions by Bosse and co-worker a'^^i'^^ and experimen- 
tal work by Imhof and Dhont'^^ for binary mixtures of 
hard spheres with very large size disparity. 

Having mentioned these antecedents, simulations 
presented here constitute, to the best of our knowl- 
edge, the first systematic investigation of slow relax- 
ation (by analyzing mean squared displacements and 
dynamic correlators) in binary mixtures covering the 
whole range between the limits of small and large dis- 
parity. They also provide a connection with results in 
Refs."^'^''^'*, which report a dynamic phase diagram dis- 
playing an A- and a B-line merging at a higher-order 
point. We suggest that the anomalous relaxation fea- 
tures observed at moderate disparity 5 ^ 2.5 for the 
system here investigated might originate from the ex- 
istence of a nearby B-line (yielding finite values for the 
non-ergodicity parameters) ending at a nearby higher- 
order point, to which anomalous relaxation features 
would be associated. The B- and A-lines would ex- 
tend from the higher-order point to, respectively, the 
limits of small [6 — I) and large size disparity. 

The article is organized as follows. In Section II 
we introduce the investigated model and give compu- 
tational details. In Section III we present simulation 
results for static structure factors, mean squared dis- 
placements and dynamic correlators. In Section IV we 
discuss dynamic quantities in the framework of the 
MCT and propose a picture for the different relax- 
ation scenarios observed by varying the size disparity. 
Conclusions are given in Section V. 

II. MODEL AND SIMULATION DETAILS 

As in the work of Ref we have simulated a mixture 
of large and small particles of equal mass m = 1 , inter- 
acting through a soft-sphere potential plus a quadratic 
term: 
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where e = 1, and a, (3 G {A, B, C, D}. Two 
sets of large particles (labelled as A and B) and of 
small particles (C and D) are introduced in order to 
avoid crystallization effects (see below). The inter- 
action is zero beyond a cutoff distance cctq/j, with 
c = 1.15. The addition of the quadratic term to the 
soft-sphere interaction, with the values Co = 7c~^^ 
and C2 = 6c~^'^, guarantees continuity of potential 
and forces at the cutoff distance. The diameters of 
the soft-sphere potential for the different types of in- 
teraction are: (Tdd 1, crcc = 1-Io-dd, o-bb = (^o-dd, 
(taa — 1-1(Tbb, and (Tq/j — (ctqq -I- app)/2 for the case 
a ^ /3. We have investigated the size disparities 6 = 
1.15, 1.6, 2.5, 5, and 8. 

The potential ^ is purely repulsive. It does not 
show local minima within the interaction range r < 
cUap- Hence, slow dynamics in the present model 
arises as a consequence of steric effects. MCT the- 
oretical works are usually carried out on systems of 
hard objects, while simulations in similar systems with 
continuous interactions are usually preferred for com- 
putational simplicity. In the present system, the tail 
of the interaction potential is progressively probed by 
decreasing temperature, which qualitatively plays the 
role of increasing packing in a system of hard spheres. 

The mixture composition is defined as the fraction 
of small particles: x^maW = {Nc + Nd)/{Na -I- A'^b + 
Nq + Nd), with Na denoting the number of particles 
of the species a. We investigate a single composition 
a^sinaii = 0.6. The number of large and small particles 
is, respectively, A^a+A^b = 1000 and Nq+Nd = 1500. 
We impose the constraints A^a = A^b and A'c = A^d- 
These constraints, together with the small selected ra- 
tios (Tcc/o'DD = ctaa/cbb = 1.1 avoid crystallization 
in the limit of small size disparity, and moreover guar- 
antee that only very small dynamic differences are in- 
duced between particles within a same set ({A,B} or 
{C,D})i. Hence, in the following we will only report 
dynamic quantities for A- and D-particles. 

Temperature T, distance, wavevector q, and time t, 
will be given, respectively, in units of e/kB, ctdd, ""ddi 
and crDD("i/e)^/^. The packing fraction is defined as: 

0=^[^AaiA + ABa|B + Aca3c-HA^D4D] (2) 

with L the side of the simulation box. We fix a value 
(j) = 0.53 for all the investigated disparities. The sys- 
tem is prepared by placing the particles randomly in 
the simulation box, with a constraint that avoids core 
overlapping. Periodic boundary conditions are imple- 
mented. Equations of motion are integrated by using 
the velocity Verlet scheme^, with a time step ranging 
from 2 X 10"'' to 5 x 10"'^, for respectively the highest 
and the lowest investigated temperature. A link-cell 
method22, is used for saving computational time in the 
determination of particles within the cutoff distance 
of a given one. 

At each state point, the system is thermalized at the 
requested temperature by periodic velocity rescaling. 
After reaching equilibrium, energy and pressure show 
no drift. Likewise, mean squared displacements and 
dynamic correlators show no aging, i.e., no time shift 
when being evaluated for progressively longer time ori- 
gins. Once the system is equilibrated, a microcanoni- 
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cal run is performed for production of configurations, 
from which static structure factors, mean squared dis- 
placements, and dynamic correlators are computed. 
For each state point, the latter quantities are aver- 
aged over typically 20-40 independent samples. 

III. SIMULATION RESULTS 

a. Static structure factors 

We compute normalized partial static structure 
factors, Sapiq) = {Pa{(l-,^)Pi3{-q,^)) I with 
Pa{q,t) — exp[zq • raj{t)], the sum extending over 
all the particles of the species a. Fig. [T]shows, for dif- 
ferent values of S, results for large-large, small-small, 
and large-small pairs. We include all the A- and B- 
particles in the set "large" and all the C- and D- 
particles in the set "small" . We present, for each value 
of S, data at the lowest investigated temperature. 

'S'largc-iargcC?) exhibits a main peak at g « 
27r((5CTDD)^^- This peak narrows and grows up by in- 
creasing S, reflecting a major ordering of the large 
particles^i^'^. Likewise, the \ow-q structure ob- 
served for small disparity dissapears. This structure 
originates from the presence of inhomogeneities or 
"voids" in the matrix of large particles, which are 
filled by the small particles. By increasing the size 
disparity, the total packing fraction is largely domi- 
nated by the contribution of the large particles, and 
the inhomogeneities progressively vanish. 

Increasing the size disparity induces opposite ef- 
fects in 5'sniaii-smaii('z)- The main peak observed at 
q « 27r CT]^]!) for 5 = 1.15 is strongly damped by increas- 
ing (5, reflecting a progressive loss of the short-ranged 
correlations for the small particles. For large dispari- 
ties, S'smaii-smaiK?) shows a nearly structureless profile 
for q > 2(7pp, close to the flat behavior expected for 
a gas. However, a peak arises at low q-values, sug- 
gesting the formation of some large-scale connected 
structure for the small particles. The presence of this 
structure should lead to an anti-correlation between 
large and small particles at the same wavevector, 
which is indeed manifested by a sharp negative peak 
in 5'iargc-smaii(<z)- Thesc structural features are some- 
what resembling of the observation of channels for 
preferential transport of alkaline ions in silica^i^i^. 
As we will show below, dynamic quantities reported 
here also present analogies with observations for the 
latter system. 

b. Mean-squared displacements 
and dynamic correlators 

Fig. [5] displays, for size disparity 6 = 1.15, the 
temperature dependence of mean-squared displace- 
ments (MSD, {[Ara{t)]'^)), density-density, and self- 
correlators for a = A- and D-particles. Density- 
density correlators for a-a pairs are calculated as 

Faa{q,t) = {pa{q,t)pa{-q,0))/{pa{q,0)pa{-q,0))- 

Self-correlators are computed as F^{q, t) — T,j exp{iq- 
[ra. j{t) — rQ,j(0)]}. The introduction of a small size 
disparity induces a small separation between the time 
scales of large and small particles. Only at the low- 
est investigated temperature, T = 0.62, this separa- 
tion is about a factor 5 in diffusivity (see Fig. 
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FIG. 1: (color online) Partial static structure factors, as 
computed from simulation data, for different size dispari- 
ties. Wavevectors are given in units of (Jj- 
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As usually observed in the proximity of liquids-glass 
transitions^, a bending occurs in the MSD after the 
initial ballistic (cx t^) regime. A plateau arises at low 
temperatures. This effect corresponds to the onset of 
the caging regime. At long times, the diffusive regime 
(cx t) is reached for values ((Ar„)^) < cr^^, i.e, when 
the particles have moved, on average, a distance of the 
order of their size. 

As usual^i^iii^, another plateau is observed for 
density-density and self-correlators (Figs. [2l3 and 
in the time interval corresponding to the caging 
regime. This interval is known as the /3-regime within 
the framework of the MCT. The correlators start to 
decay from the plateau at times corresponding to the 
onset of the diffusive regime in the MSD. This sec- 
ond decay is known as the a-regime, and it is often 
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FIG. 2: (color online) Symbols: simulation results for size 
disparity 5 — 1.15 at different temperatures. Two identi- 
cal symbols (empty and filled for respectively A- and D- 
particles) correspond to a same temperature. Panel (a): 
mean-squared displacements; panel (b): density-density 
correlators; panel (c) : self-correlators. The two latters are 
plotted at a fixed wavevector q — 5.56a^^. Straight lines 
in panel (a) correspond to linear behavior. 



described by an empirical Kohlrausch- Williams- Watt 
(KWW) function, given by a stretched exponential, 
Aqexp[—{t/Tq)l^'>], with Aq < 1 the plateau height 
and < /3g < 1. The parameters Aq, Tq and Pq are 
g-dependent. 

Fig. [3] displays results for disparity 5 = 1.6, which 
induces a clear separation between the time scales for 
A- and D-particles. Interestingly, the MSD for D- 
particles shows, at sufficiently low temperatures, an 
apparent sublinear power-law behavior, oc f", over 



FIG. 3: (color online) As Fig. [2]for 5 — 1.6. The wavevec- 
tor for panels (b) and (c) is g = 4.520-^^. Straight lines 
in panel (a) correspond to linear or sublinear power-law 
behavior (exponents are given). 



time intervals of almost two decades after the caging 
regime. Hence, differently for the standard behavior, 
the small particles need to move distances of a few 
times their size to reach the diffusive regime. The 
exponent a < 1 seems to decrease by decreasing tem- 
perature. All these features suggest that increasing 
the size disparity induces a progressive separation in 
the time scales of the large and small particles, and 
a qualitative change in the relaxation scenario for the 
small particles. 

These trends are confirmed by results for S = 2.5 
(see Fig. |4|). D-particles reach the diffusive regime 
only for displacements of about ten times their size. 
Density-density correlators for the D-particles exhibit 
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FIG. 4: (color online) As Fig. [3]for S = 2.5. The wavevec- 
tor for panels (b) and (c) is g = 2.45(7^]^. The straight 
line in panel (b) indicates logarithmic-like behavior. 



extreme stretching. No defined plateau is present. 
At the lowest investigated temperatures a logarith- 
mic decay is observed over a time window of about 
four decades. Self-correlators for the D-particles show 
a qualitatively different behavior. They exhibit a de- 
cay much faster than density-density correlators at the 
same temperature. This feature is displayed in Fig. [5] 
by comparing, at a fixed temperature T = 0.30 and 
several wavevectors, results for Fdd(9, t) and t). 
Both correlators only converge in the limit of high- 
q. The time scale separation at moderate and low 
wavevectors reaches values of even two decades. It 
is interesting to note that, after a fast decay down 
to Fl,{q,t) ~ 0.3, self-decorrelation is considerably 
slowed down, as evidenced by the long tail observed at 



FIG. 5: (color online) Dynamic correlators of D-particles, 
at different wavevectors, for size disparity 5 = 2.5 and T = 
0.30. Empty and filled symbols correspond, respectively, 
to self- and density-density correlators. 



long times (Fig. Hj:). This tail is observed in the same 
time scale where the MSD exhibits power-law subd- 
iffusive behavior, suggesting a common microscopic 
origin for both features. 

A new relaxation scenario arises by increasing the 
size disparity to much larger values. Fig. [6] displays 
simulation results for 6 = 5. In order to provide a 
clearer visualization of the decay for density-density 
and self-correlators, data for the latters are plotted as 
a function of tT^^^, i.e., rescaling the time by the ther- 
mal velocity. As expected, the large particles again 
exhibit a much slower dynamics. Density-density cor- 
relators for the small particles do not show a slow 
decay at any temperature, even for wavevectors corre- 
sponding to wavelengths of several particle diameters. 
Decorrelation occurs in an essentially exponential way 
down to very small values of i^DD {q, t) , where a nearly 
flat background arises below some temperature and 
extends to very long times. This feature is progres- 
sively enhanced by increasing 5, as shown in Fig. [7] 
for (5 = 8. The amplitude of the background increases 
as the system is cooled down. Qualitatively similar 
results are observed for other wavevectors. 

Self-correlators for the small particles display a 
striking result. As for density-density correlators, 
they exhibit a fast decay followed by a long tail. How- 
ever, they decay to zero at all the investigated temper- 
atures (Figs. [H; and[71;), even at those where i^DD(9, t) 
shows a finite long-time limit. Results for the MSD in 
Figs. [5^ and[7^ evidence that the small particles reach 
the long-time diffusive regime at such temperatures. 
All these features demostrate that the self-motion of 
small particles is ergodic at all the investigated tem- 
peratures. A small particle can reach regions arbi- 
trarely far from its initial position. However, accord- 
ing to results in Figs, ^jp and [7|d, coherent dynamics 
are non ergodic below some given temperature. 

This unusual feature and its precursor effects shown 
above for 6 = 2.5 can be tentatively assigned to the 
existence of preferential paths for the motion of the 
small particles, forming a slowly relaxing large-scale 
structure. As a consequence, small particles would be 
able to diffuse along such a structure, but the associ- 
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FIG. 6: (color online) As Fig. [3]for 5 = 5. The wavevector 
for panels (b) and (c) is g = 2.06(7]^^. Time in panels (b) 
and (c) is rescaled by T^^^. The horizontal dashed lines 
indicate the zero value. 



ated density fluctuations would be long-lived, leading 
to very slow coherent dynamics as compared to self- 
motions. At very large size disparities the matrix of 
large particles becomes glassy and the paths for dif- 
fusion of the small particles would not be able to re- 
lax, leading to non-ergodic coherent dynamics while 
self-motions remaining ergodic. Though the charac- 
terization of this eventual structure of preferential 
paths for the small particles is beyond the scope of 
this article, the presence of low-g peaks in the par- 
tial static structure factors (Figs. [T]d and [T}:) sup- 
ports this hypothesis. It is worth mentioning that a 
time scale-separation for incoherent and coherent dy- 
namics similar to that reported here for S = 2.5 has 



FIG. 7; (color online) As Fig. |6]for 5 = 8. The wavevector 
for panels (b) and (c) is q = 0.65a^^. 

been observed for the diffusion of alkaline ions in sil- 
ica matrixes^. In this system a channel-like structure 
for preferential diffusion of the alkaline ions has been 
explicitely characterized^^ji^iii. The separation of in- 
coherent and coherent dynamics for the formers has 
been assigned to the exsitence of this structure. 

Results reported in this section evidence that in- 
creasing the size disparity yields unusual relaxation 
scenarios for the dynamics of the small particles. In 
the next section we discuss these results within the 
framework of the MCT and propose a global physical 
picture for the different scenarios observed at different 
disparities. 

IV. DISCUSSION 
a. MCT framework 
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Next we summarize the basic predictions of MCT 
and test them in the present system. In its ideal ver- 
sion, which neglects activated hopping events, MCT 
predicts a sharp transition from an ergodic liquid to 
a non-ergodic arrested state (glass) at a given value 
of the relevant control parameter (in the following the 
temperature, 7^) ^,8,10 ^ When crossing the transition 
point T = Tc (critical temperature) from the ergodic 
to the arrested state, the long-time limit of the corre- 
lator for wavevector q, F{q,t), jumps from zero to a 
non-zero value, denoted as the critical non-ergodicity 
parameter, Moving beyond the transition point 
into the non-ergodic state yields a progressive increase 
of the non-crgodicity parameter, fq > In the 

standard case the jump in F{q,t) is discontinuous, 
i.e., the critical non-ergodicity parameter takes a 
finite value. MCT transitions with > are denoted 
as type-B transitions. For ergodic states close to the 
transition point, the initial part of the a-process — 
i.e., the von Schweidler regime — is given by a power 
law expansioniiSiiiS: 



2b 



(3) 



(2) 

with < 6 < 1. The prefactors hq and hq only de- 
pend on the wavevector q and are different for each 
correlator. The characteristic time Tc only depends 
on the separation parameter \T ~Tc\ and is divergent 
at the transition point. As mentioned above, the de- 
cay from the plateau to zero is often described by an 
empirical KWW function with a q-dependent stretch- 
ing exponent /3q. An interesting prediction of MCT— 
is that Pq = b in the limit of large-g. This result 
has been widely tested^'^1^'^1^'^'^'^ and provides 
a consistency test for data analysis. 

Another prediction of MCT for state points close 
to the transition point is the power law dependence of 
the diffusivity and the a-relaxation time: 
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The a-relaxation time, Tq,, is a time scale probing the 
a-process. In practice it can be defined as the time 
where F(q, t) decays to some small value a:, provided 
it is well below the plateau. The exponent 7 in Eq. 
is given hy^j^^: 
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with < a < 0.4. Hence 7 > 1.75. The critical 

exponents a, 6, and 7 are univoquely related with the 
so-called exponent parameter 1/2 < A < 1 through: 



A = 



r^(i + b) _ r^(i~a) 

r{l + 2b) " r(l -2a)' 



(6) 



FIG. 8: (color online) MCT analysis for 5 = 1.15. Sym- 
bols are simulation results. Panel (a): For T = 0.67, q- 
dependence of the stretching exponents obtained from fit- 
ting the decay from the plateau to a KWW function. The 
corresponding correlators are given in the legend. The 
dashed line indicates the limit I3{q —> 00) ~ 0.48. Panel 
(b): density-density correlators for D-D pairs, at T = 0.67, 
for different wavevectors. Curves are fits to Eq. with 
b — 0.48. Panel (c): diffusivities and inverse relaxation 
times for q = Q.Za^^ (see text for the definition). Empty 
and filled symbols correspond, respectively, to A- and D- 
particles. Dashed and solid straight lines are fits to Eq. Q 
for, respectively. A- and D-particles. The values Tc = 0.64 
and 7 = 2.8 are forced. 



with F the Gamma functioniiSiiiS. The exponent pa- 
rameter A is univoquely determined by the static cor- 
relations (i.e., by the total and partial static structure 
factors) at the transition point T = T^. 

When numerical solutions of the MCT equations 
are not available the non-ergodicity parameters, pref- 
actors and exponents in Eqs. pi4l5l6p are obtained as 



fit parameters from simulation or experimental data. 
Consistency of the data analysis requires that the so- 
obtained values for the exponents fulfill Eqs. ([5]) and 
([6|). It must be stressed that values of MCT critical 
exponents obtained as fit parameters must be taken 
with care, since the range of validity (in time and tem- 
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FIG. 9: (color online) As Fig. [8]: for 5 = 2.5. The values 
Tc — 0.19 and 7 — 4.0 are forced. The wavevector for the 
Qf-relaxation times is q — 2.4:5a^ly. 



perature windows) of the asymptotic Eqs. (|3|4|) is not 
known a priori. The case of the power law (j4]) is spe- 
cially problematic, not only for the upper limit in tem- 
perature, but also for the lower one. Deviations of the 
power-law behavior are often observed below some ill- 
defined temperature very close to ^j^g 
presence of activated hopping events which are not ac- 
counted by ideal MCT. Hence, by selecting different 
temperature ranges (both the upper and lower lim- 
its being ill-defined) for fits to Eq. ^ one can find 
different effective exponents 7, which can be rather 
different from the actual one. It often occurs that the 
effective exponents obtained for the diffusivities are 
very different from those obtained for the relaxation 
times^i^i^. See, e.g., Refsi^!^ for a discussion on 
this point. 

Fig. [8] shows a test of MCT predictions for size 
disparity S = 1.15. Within numerical uncertainties, 
the stretching exponents of self- and density-density 
correlators, both for large and small particles, seem to 
approach a common value f3 « 0.48 in the large-g limit 
(Fig. [5]a,), in agreement with MCT expectations^. 
Consequently with the prediction P{q — > 00) = b, we 
fix the von Schweidler exponent b to this value and fit, 
for the different correlators, the initial decay from the 
plateau to Eq. Consistently, a good description 

of the simulation data is achieved. Fig. [SJd shows the 
corresponding fits for FddC'ZjO- 

The value b — 0.48 provides through Eqs. ([5]) and 
dSl) the values A = 0.79, a = 0.28, and 7 = 2.8. 
Now we test the validity of Eq. (|4]) with this latter 
value of 7. Fig. ^ shows fits to Eq. ^ of diffu- 
sivities and Qf-relaxation times. The diffusivities Da 
are calculated as the long-time limit of {[Ara{t)]'^) /6t. 
We use X = 0.2 for the definition of the a-relaxation 
times (denoted as for self-correlators and Taa for 
density-density correlators). The latters are shown for 
q — 6.3(Tpj!). We have fixed b = 0.48 and forced a com- 
mon Tc for all the diffusivities and relaxation times, 
as predicted by MCT for monodisperse spheres or bi- 
nary mixtures with small size disparity^. The best 
fits are achieved with Tc = 0.64. As expected, the 
power law ^ fails above some high temperature. For 
the diffusivities it also fails at temperatures very close 



FIG. 10: (color online) Density-density correlator for A-A 
pairs at disparity S = 2.5 and temperature T — 0.60. Dif- 
ferent symbols correspond to different wavevectors. The 
straight line indicates logarithmic decay. 



to Tc, where activated hopping events are expected to 
raise the simulation values above the theoretical pre- 
dictions. It is often found that these deviations arise 
for diffusivities at higher temperatures than for relax- 
ation times^. This is not surprising, since diffusivities 
and relaxation times are dominated by the contribu- 
tion of, respectively, fast and slow particles. Similarly, 
deviations due to hopping events in binary mixtures 
are stronger for the small particles^i*^'^ , which are, 
on average, faster than the large ones. This feature 
is indeed observed in Fig. [5J: and is enhanced by in- 
creasing 6 (see Fig. |9]for 5 = 2.5), suggesting a major 
relevance of hopping events for the diffusivity of the 
small particles at larger disparity. As discussed above, 
a fully free fit to Eq. Q would provide rather differ- 
ent exponents 7 for relaxation times and diffusivities. 
Only a simultaneous analysis of dynamic correlators, 
as illustrated in Figs. [S^ andjHb can provide a robust 
determination of the critical exponents. 

An analogous analysis for disparity S = 1.6 pro- 
vides for the critical exponents the values a = 0.25, 
b = 0.40, 7 = 3.2, and A = 0.85. We obtain a crit- 
ical temperature Tc = 0.42, much lower than that 
for 6 — 1.15. Hence, increasing the size disparity 
to S = 1.6 at constant packing fraction shifts the 
glass transition to lower temperatures, as predicted 
by MCT^^ for similar values of S and observed in 
experiments^ and simulations^!^. 

This trend is also fulfilled for S — 2.5, for which the 
critical temperature is Tc — 0.19. Fig. [H] shows sim- 
ulation results for the diffusivities and the relaxation 
times. Due to the small amplitude of the tail exhib- 
ited by the self-correlators for the D-particles (see Fig. 
[It) , we have defined the corresponding relaxation time 
Tx by using x = 0.03. The value x = 0.2 is used for 
the other relaxation times displayed in Fig. [9l We 
obtain the values a = 0.21, b = 0.31, 7 = 4.0, and 
A = 0.90 for the critical exponents^S. Another pre- 
dicted trendS that is fulfilled by increasing S is the 
observed increase of the exponent parameter A. 

In our previous article (see Fig. 9 in Ref,— ) we have 
analyzed the decay of FDD(g, t) a.t 6 = 2.5 in terms of 
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a logarithmic expansion, 

F{q, t) « ~ Hg ln(i/Tc) + i/f ) In2(t/Te), (7) 

instead of the von Schweidler law ([3]). Within the 
framework of MCT, Eq. ([7]) is associated to a nearby 
higher-order transitio n^^i^^'^^ , which is characterized 
by a value of the exponent parameter A = 1, or to a 
standard transition with a value of A very close to 1. 
Hence, the value A = 0.90 here obtained for disparity 
5 = 2.5 justifies the description used in Ref.— in terms 
of a logarithmic expansion. 

It must be stressed that, in principle, this choice 
is not in contradiction with the description, for this 
same value of S — 2.5, of the correlator for the A- 
particles in terms of a von Schweidler law, as done 
in Ref.'' (see Fig. 7b there). According to MCT the 
presence of a nearby higher-order transition is mani- 
fested in a given correlator F(q, t) by a purely loga- 
rithmic decay at some given temperature. Moreover, 
by varying q at that temperature, the shape of the 
decay exhibits a convex-to-concave crossove r^^i^^i^°i^^ . 
We have observed these features for the small parti- 
cles at low temperatures (see Fig. 9a in Ref4), but 
they are absent for the large particles at the same 
temperatures (Fig. 7b in Ref."'). However, they are 
clearly observed at higher temperatures, as shown 
here in Fig. [10] for T — 0.60. Logarithmic relax- 
ation covers more than two time decades at q w 5(7^^. 
The fact that these anomalous features (logarithmic 
relaxation and convex-to-concave crossover) are ob- 
served for the large and the small particles at dif- 
ferent temperatures is not, in principle, a failure of 
MCT. The optimal distance to the transition point for 
the observation of anomalous relaxation is determined 
by the coefficient Hq"^^ in Eq. ([7]). This coefficient 

is decomposedi^ as = ^({x"}) + B({x"})7<'q, 

where the g-independent terms A({x"}) and i?({x"}) 
depend on the state point {x"} = {T, (5, Xsmaii, 0, ■■•} 
in the control parameter space, and Kg, which is de- 
termined by static correlations, only depends on q and 
is different for each correlator. As a result of this de- 
composition, for a given correlator there are paths in 

(2) 

the control parameter space where Hq ~ and along 
which, according to Eq. ([7]), a purely logarithmic de- 
cay will be observed. The fact that in the present 
case the partial static structure factors for the large 
and the small particles are very different (Fig. [T]) 
might yield very different i^Tg-coefficients for the re- 
spective correlators, and as a consequence, very dif- 
ferent paths for purely logarithmic relaxation of the 
large and the small particles. This feature might ex- 
plain why anomalous relaxation for both types of par- 
tices is observed at very different temperatures. Un- 
less one follows the optimal path in the control pa- 
rameter space — which in principle involves a simul- 
taneous variation of several control parametersi^iiS — 
anomalous relaxation vanishes by decreasing temper- 
ature, and a standard two-step decay is recovered (see 
Refi^ for an illustrative example). Hence, correlators 
for the large particles at low temperature are well de- 
scribed by a von Schweidler law. The fact that anoma- 
lous relaxation is present for the small particles at very 
low temperatures suggests that decreasing tempera- 
ture and fixing the other control parameters does not 
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FIG. 11: (color online) MCT analysis for 3^5. Sym- 
bols are simulation results. Panel (a): For T = 0.50, 
g-dependence of the stretching exponents obtained from 
fitting the decay from the plateau to a KWW function. 
The corresponding correlators are given in the legend. The 
dashed line indicates the limit P{q —> oo) ~ 0.54. Panel 
(b): density-density correlators for A-A pairs, at T = 0.50, 
for different wavevectors. Curves are fits to Eq. (|3} with 
b = 0.54. Panel (c): diffusivities and inverse relaxation 
times for q — 1.33a^l^. Empty and filled symbols corre- 
spond, respectively, to A- and D-particles. Dashed lines 
are fits to Eq. JU. The values Tc = 0.42 and 7 = 2.6 are 
forced. 



deviate the system far from the corresponding optimal 
path for the small particles. 

The analysis for the large particles ai S — 5 (see 
Fig. [TT]) provides a critical temperature Tc = 0.42. 
Hence, the trend observed for smaller disparities is re- 
versed for sufficiently large values of S, and the glass 
transition exhibits a reentrant behavior by using 6 as 
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FIG. 12: (color online) Circles: Critical temperatures Tc 
for the different investigated disparities. Squares are the 
corresponding values of the exponent parameter A. Note 
the logarithmic scale in 5. Error bars are given by the 
symbol size. 



a control parameter. Increasing disparity at constant 
packing fraction initally decreases Tc. However, the 
latter reaches some minimal value and beyond some S 
it starts to increase. For 6 — 5 we obtain the critical 
exponents a = 0.30, b = 0.54, 7 = 2.6, and A — 0.76. 
Hence, the reentrant behavior observed for Tc also in- 
volves a non-monotonic behavior for the exponent pa- 
rameter A. Increasing d initially raises A to large val- 
ues (0.90 for S = 2.5), suggesting an eventual nearby 
higher-order transition (A = 1). Beyond some dis- 
parity, for which it reaches a maximum, the exponent 
parameter starts to decrease. These trends are con- 
firmed for the largest investigated disparity 6 = 8, for 
which we obtain Tc = 0.52 and the critical exponents 
a = 0.31, b = 0.57, 7 = 2.5, and A = 0.74. We note 
on passing that data for correlators of the large parti- 
cles in Figs, ini and [7| show a decay from the plateau at 
temperatures far below Tc. This feature is in principle 
related to the mentioned hopping events not included 
in the ideal version of MCT, which restore ergodic- 
ity below Tc. The fact that, despite the final decay, 
the plateau height exhibits a clear increase below the 
value estimated for Tc is consistent with expectations 
of ideal MCT. 

Fig. [T2] displays the values of Tc and A obtained 
at the different investigated disparities. Only small 
variations of these values provide reasonable descrip- 
tions of simulation results in a consistent way (i.e., by 
fitting simultaneously different data sets as exposed, 
e.g., in Fig. [5]). The so-estimated error bars for Tc 
and A are given by the symbol size in Fig. [121 The 
values of A obtained for 6 < 1.6 are somewhat higher 
than those predicted by MCT^^ in the same range for 
a binary mixture of hard spheres at the same com- 
position Xsmaii = 0.6. From Fig. 13 in Ref.^^ (note 
the transformations of Xsmaii and d to the quantities 
represented there), we estimate A w 0.73 for S ~ 1.15 
and A « 0.78 for S = 1.6. Unfortunately, values of A 
for larger disparities have not been reported. 

It must be stressed that determining wether there 
is actually a higher-order point (A = 1) at some dis- 
parity around 5 = 2.5, or wether A instead reaches a 



maximum 0.90 < Xmax < 1, is a difficult task. Since, 
from Eq. ([S]), A — *■ 1 implies that 6 ^ 0, a deter- 
mination of the specific value of b through the limit 
b = (3{q — > 00) is not reliable for extremely small val- 
ues of b. Only the solution of the MCT equations can 
unambiguously determine the actual value of Amax- 
Still, the existence or absence of a higher-order point 
is not a question so important as the anomalous re- 
laxation features associated to it (logarithmic decays 
and convex-to-concave crossover in dynamic correla- 
tors, as well as sublinear power-law behavior for the 
MSD^^), which are also expected for values A 1 as 
those found here. 

For the case of the small particles at disparities 
5 = 5 and 5 = 8 we have not performed an analysis in 
terms of MCT similar to those above exposed. Indeed, 
data for the small particles in Figs. [S] and [7] exhibit 
a very different scenario of that observed for 6 < 2.5. 
As mentioned in Section HI, density-density correla- 
tors do not show a slow decay, but a fast one followed 
by a nearly flat background below some temperature, 
its amplitude growing up by decreasing temperature. 
We identify the background amplitude as an opera- 
tional non-ergodicity parameter that increases from 
some critical value As shown in Figs. [HId and[7j5 it 
is difficult to establish wether defined in this way is 
zero, but it is clear that it takes, as much, a extremely 
small value. We can also define an operational Tc for 
the small particles as the temperature where starts 
to increase. Though a precise determination of this 
temperature is difficult, simulation data suggest that 
0.4 < Tc < 0.75 and 0.5 < Tc < 0.75 for, respec- 
tively, (5 = 5 and 6 = 8, i.e., of the order of the critical 
temperature for the large particles (see above), and 
compatible with a common Tc for all the particles. 

Data for 6 = 5 show a final decay of the background 
to zero (see Fig. |BJd). This decay is, again, presumably 
related with the presence of hopping events. Note that 
despite of this decay, the increase of the background 
amplitude for T < 0.4 is evident, as expected for the 
crossing of a MCT critical temperature. 

The behavior reported here for density-density cor- 
relators of the small particles at 5 = 5 and 5 = 8 re- 
sembles features characteristic of MCT transitions of 
the so-called type-A, which are defined by a zero value 
of the critical non-ergodicity parameter, in contrast 
to the finite value defining the standard type-B tran- 
sitions. A recent realization of this scenario has been 
reported for a system of dumbbell molecules^'^'^. 
While for moderate molecular elongations a standard 
relaxation scenario is observed, for small elongations 
angular correlators of odd order exhibit features anal- 
ogous to those of Figs. [6)3 and [7):;^^. Theoretical cal- 
culations for this latter system^^ relate such features 
to the existence of a MCT transition of the type-A. 

Data presented here are consistent with early MCT 
theoretical calculations by Bosse and co-worker a'^^i'^^ 
and later by Harbola and Das^^ for a binary mix- 
ture of hard spheres with the packing fraction as the 
control parameter. Though in these works no in- 
formation is given about relaxation features, calcu- 
lations are reported for non-ergodicity parameters of 
self- (/p and density-density (fq) correlators. These 
calculations predict, for sufficiently large disparities, a 
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glassy phase for density-density correlations (/, > 0) 
where self-correlations for the small particles remain 
ergodic (/| = 0)'^^ and the diffusivity for the latters 
is finiteS^, in agreement with simulation results pre- 
sented in this article and with early experiments by 
Imhof and Dhont^ in binary mixtures of colloidal 
silica particles with size disparity S = 9.3. While 
density-density correlators of the small particles dis- 
play a critical temperature Tc > 0.4, Figs. [SJ: and 
[71; do not show, even for T <^ T^^ any signature of a 
crossing of a critical temperature for self-correlators. 
Likewise, small particles exhibit large values of the dif- 
fusivities (Z?d > 0.02 for (5 = 8) and very weak caging 
effects in the MSD (Figs. [6^ and[7K) at temperatures 

As mentioned in Section III, transport of alkaline 
ions in silica display a separation between the time 
scales for coherent and incoherent dynamics similar 
to that presented here for 5 = 2.5. In a recent work, 
MCT calculations by Voigtmann and Horbach have 
reproduced this feature^'*, which has also been in- 
terpreted as a precursor effect of an eventual type- 
A transition^. Given the similarities (already noted 
in^) of this system with binary mixtures of suffi- 
ciently large size disparity, the former might be a can- 
didate to display, by properly tuning the concentra- 
tion of alkaline ions, many of the unusual relaxation 
features reported here. Another possible candidate is 
a mixture of star polymers with arms of very different 
length. In a coarse-grained picture, this system can 
be seen as binary mixture of ultrasoft sphere o'^'^^ with 
large size disparity. By properly tuning the mixture 
composition and the disparity it is possible to obtain 
a glassy matrix formed by the large stars where the 
small ones remain fluid'^. 

To finish this subsection, we make a brief comment 
on MCT predictions for Dt, the product of the dif- 
fusivities by the a-relaxation time. Fig. [T2] shows, 
for different size disparities, results for Dt as com- 
puted from numerical data of Figs. [8t, [H and [TTb . 
This product approaches a constant value at high tem- 
peratures and, for some of the data sets in Fig. [T51 
increases by orders of magnitude by decreasing tem- 
perature. This feature has been widely observed in 
very different systemsi^i^i^i^i^i^iS and is usually 
attributed to the presence of dynamic heterogeneities, 
which enhance the diffusivity (dominated by the con- 
tribution of the fast particles) as compared to the a- 
relaxation time (dominated by the contribution of the 
slow particles). It must be emphasized that devia- 
tions of Dt from constant behavior are not a priori a 
signature of deviations from MCT predictions. Since 
the power law D cc {T — T^)'^ is just a leading-order 
result, MCT does not predict a temperature indepen- 
dent behavior for the product Dt, but just a weak 
variation between the high-T limit and a bounded 
value at T ^ Tc (see e.g., Refs,^iS). This predic- 
tion is compatible with data for the large particles at 
5 = 2.5 (Fig. [T3d) and (5 = 5 (Fig. [H;). Indeed 
only small deviations from MCT, in the investigated 
T-range, are observed by looking separately at Z?a, 
taa, and ta (Figs. [9l and [TTt). Deviations for the 
diffusivity are larger at 5 = 1.15 (Fig. [5t), leading to 
an increase of about an order of magnitude in DaJaa 
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FIG. 13: (color online) For different size disparities, prod- 
uct of diffusivities D by a-relaxation times r. The values 
of the latter quantities, as well as the critical temperatures 
Tc, are the same ones as in Figs. |8};, [H and lllb . All the 
axes are in logarithmic scale, except for the vertical one in 
panel (c), which is in linear scale. 



and Data within the investigated T-range (Fig. [T3k). 
On the contrary, Dt for the small particles strongly 
deviates from the MCT prediction in all cases, as ex- 
pected from the strong deviations already observed by 
looking separately at D^. 

b. A global picture 

The different relaxation scenarios here reported for 
the small particles by varying the size disparity pro- 
vide a connection with MCT theoretical results by 
Krakoviack in Refs.— for a mixture of fixed and 
mobile hard spheres, a model introduced to reproduce 
qualitative dynamic features of fluids confined in ma- 
trixes with interconnected voids. In that work, the 
dynamic phase diagram of the mobile particles dis- 
plays two lines of type-A and type-B transitions in the 
plane 0fix-0mob, where (pax and (/'mob are the packing 
fractions of, respectively, the fixed and mobile parti- 
cles. The A- and B-lines merge at a higher-order point 
(namely an Aa-poinl;^). The B-line extends from the 
Aa-point to the limit (/)fix — 0, where the fluid of hard 
spheres is recovered. The A-line extends from the A3- 
point to the Lorentz gas limit (a single particle diffus- 
ing in a matrix of fixed obstacles) at (pmoh = 0. 

For high concentrations of the mobile particles, the 
matrix of fixed particles does not yield signifhcant con- 
finement effects. The only caging mechanism is nor- 
mal hard-sphere repulsion at short length and time 
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scales, and dynamic correlators exhibit a standard 
two-step decay^. The transition point is of the type- 
B and hence the jump of the long-time limit of the 
density-density correlator is finite^*, providing a non- 
zero value of the critical non-ergodicity parameter. 

For high dilution of the mobile particles, hard- 
sphere repulsion at short length and time scales 
does not yield temporary caging. As a consequence, 
density-density correlators display a fast decay to val- 
ues close to zero^i. At longer time scales the mobile 
particles probe the structure of the confining matrix of 
fixed particles, which leads to a "mesoscopic" caging, 
characterized by a length scale larger than that char- 
acteristic of bulk- like hard-sphere repulsion. As a con- 
sequence of this large-scale caging mechanism (con- 
finement), density-density correlators for the mobile 
particles exhibit a long tail of small amplitude after 
the fast microscopic decay^i. At the transition point, 
of the type-A, the long-time limit does not exhibit a 
finite jump but grows up continuously'^'*, providing a 
zero value for the critical non-ergodicity parameter. 

As mentioned above, at moderate concentrations 
of fixed and mobile particles, a higher-order Aa-point 
arises^ as a consequence of the competition between 
the mentioned "microscopic" and "mesoscopic" caging 
mechanisms. Relaxation features have not been re- 
ported in Refsi^^i^i for state points close to the A3- 
point. However, as stressed in^^, they will necessarily 
display the anomalies reported here (logarithmic de- 
cay and convex-to-concave crossover in dynamic cor- 
relators, and sublinear power-law behavior for the 
MSD), as a mathematical consequence of the value 
A = 1 defining the Aa-point^. 

How do results in Refsi^i^i compare with relaxation 
features presented here?. First it is worth emphasiz- 
ing that, differently from the mobile particles in that 
work, high dilution of the small particles is not the 
key ingredient for yielding a type-A transition for the 
latters in the system here investigated. Indeed, data 
reported in our previuos article^ for S = 2.5, (p = 0.53, 
and Xsmaii = 0.1 show features rather different from 
those characterizing type-A transitions. It must be 
noted that, for these control parameters, the system 
here investigated is much denser than at state points 
close to the A-line for the mixture of fixed and mobile 
particles of Refi^. Despite its low density, decorre- 
lation of the mobile particles in the latter system is 
blocked at large length scales due to the absence of 
percolating free volume. In the system here investi- 
gated, relaxation is possible at higher densities due 
to the non-static nature of the confining matrix. The 
slow motion of the large particles creates regions of 
sufficient local free volume which facilitate decorrela- 
tion of the small particles. As a consequence of high 
density, short-range bulk-like caging is a relevant ar- 
rest mechanism for small particles and leads to a slow 
decay of dynamic correlators (Fig. 10a in Ref.*), dif- 
ferent from features characteristic of type-A transi- 
tions. Hence, for not sufficiently large disparities as 
d — 2.5, type-A transitions cannot exist at any mix- 
ture composition. 

A way of yielding a type-A scenario for the small 
particles in the present system — at a fixed pack- 
ing fraction and mixture composition — is by increas- 




FIG. 14; (color online) Typical slab (of thickness IOctdd) 
at T = 0.03 for size disparity 5 = 8. Dark: A- and B- 
particles. Light: C- and D-particles. 
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FIG. 15: (color online) Symbols: non-ergodicity parame- 
ters of Fou {q, i) for different size disparities. Results for 
5 = 1.15 and 5 = 2.5 are the critical values. Results for 
5 = 8 correspond to a state point, T = 0.03, below the A- 
line (see text). Note the double logarithmic scale. Dashed 
lines are guides for the eyes. 



ing the size disparity, as shown above for 5 — b and 
(5 = 8. In this situation small particles move in a 
medium of low local density and short-range bulk-like 
caging is supressed. This effect is illustrated in Fig. 
[HI which depicts a typical slab of the simulation box 
for (5 = 8. However, small particles are caged at large 
length scales where they probe the structure of the 
confining medium. This crossover from "microscopic" 
to "mesoscopic" caging for the small particles by in- 
creasing the size disparity is manifested by a strong 
narrowing of the g-dependence of the non-ergodicity 
parameter (see Fig. I15p . which is a signature of an 
increasing localization length'^^ . Below a certain tem- 
perature, where the tail of the soft-sphere potential 
is probed, the effective packing of the confining ma- 
trix will become sufficiently large to freeze density- 
density correlations for the small particles, leading to 
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a dynamic picture analogous to the type-A scenario 
observed for the system of RefsjSiiS^. 

Relaxation features reported here for the small par- 
ticles at 5 = 5 and 5 = 8 suggest the presence of 
a nearby A-line, originating from large-scale caging 
(confinement) induced by the slow matrix of large 
particles. For small size disparity, S 1, one obvi- 
ously recovers the standard MCT scenario predicted 
for the monoatomic fluid of hard spheres and widely 
observed in mixtures of hard^, soft^S, or Lennard- 
Jones spheres^, where small disparity is just intro- 
duced as a way of avoiding crystallization. The stan- 
dard scenario, originating from short-range bulk-like 
caging, is characterized by a nearby B-line. 

Increasing the size disparity at a fixed mixture com- 
position will weaken the effects of bulk-like caging and 
strengthen those associated to confinement. This fea- 
ture suggests that a crossover from a B- to an A-line 
will occur by increasing the value of (5, in analogy with 
results for the mixture of fixed and mobile spheres 
investigated in Refsi^s^i, where such a crossover is 
obtained by varying the composition. For the latter 
system, the B- and A-lines merge at a higher-order 
point. The resemblance of results reported here and 
in Ref.— for 5 = 2.5 with relaxation features character- 
izing higher-order MCT transitions suggest an analo- 
gous merging. The observed variation of the exponent 
parameter A by increasing 6 (Fig. [T2)) also supports 
this hypothesis, or at least, the existence of a "quasi" 
higher-order point with A very close to 1. Anoma- 
lous relaxation at 6 ^ 2.5 would originate from the 
presence of a nearby B-line (due to the finite value 
of the critical non-ergodicity parameters), ending at 
a nearby higher- (or quasi-higher) order point, which 
would produce the observed anomalous relaxation fea- 
tures. 

V. CONCLUSIONS 

We have carried out simulations on a mixture 
of large and small particles. Slow dynamics have 
been investigated for a broad range of size dispar- 
ities at a fixed mixture composition, and discussed 
within the framework of the Mode Coupling Theory 
(MCT). Results presented here extend a recent in- 
vestigation on the composition dependence at a fixed 
large disparity^. Increasing the size disparity yields a 
crossover between different relaxation scenarios. For 
sufficiently large disparities an anomalous relaxation 
scenario arises, displaying features very different from 
the standard ones observed at small disparity. Some of 
these features, which resemble predictions of MCT for 
state points close to higher-order transitions, are sub- 
linear power-law behavior for mean squared displace- 
ments, and logarithmic decays and convex-to-concave 
crossovers for dynamic correlators. 

For very large disparities, a new scenario arises, 
showing features associated to MCT transitions of the 



type-A, which are characterized by a zero value of 
the critical non-ergodicity parameter, different from 
the finite value defining the B-transitions observed at 
smaller disparities. In this scenario, self-correlations 
for the small particles remain ergodic at tempera- 
tures far below the freezing of density-density correla- 
tions. Small particles remain fluid in the glassy matrix 
formed by the large particles. 

All these features provide a connection with MCT 
theoretical results for a mixture of mobile and fixed 
particles^^^, which report a dynamic phase diagram 
displaying an A- and a B-line merging at a higher- 
order point. A similar crossover is suggested for the 
system here investigated by varying the size disparity. 
If this hypothesis is correct, the global MCT picture 
discussed here might not be a specific feature of liquids 
confined in matrixes with interconnected voids^^i^, 
but a more general one. Of course, a proper answer 
to this question can only be provided by solving the 
corresponding MCT equations. However, the highly 
non-trivial observed analogies suggest to consider it 
as a plausible hypothesis and might motivate theoret- 
ical work in this problem within the framework of the 
MCT. To our knowledge, there is only an approach 
different from MCT which has been able to reproduce 
the basic phenomenology characterizing the "higher 
order-like" scenario. Namely a dynamic facilitation 
picture, which has reproduced reentrant behavior of 
diffusivities for short-ranged attractive coUoids^^ as 
well as concave-to-convex crossovers^ and logarith- 
mic decays in dynamic correlator o^^'^^ , and sublinear 
power-law behavior in the MSE>S. However, a simi- 
lar picture is missing for the "higher order-like" sce- 
nario as an intermediate state of a crossover between a 
"type B-like" and a "type A-like" scenario by varying 
the size disparity. 

The mixture of soft spheres investigated in this 
work shares many dynamic features with other mixed 
systems of different nature in the case that the two 
components show very different mobilities, as poly- 
mer blendsSi, colloidal mixtures^i^, star polymer 
mixtures^, or ion conducting glasses^^*^. Such fea- 
tures in this latter systems have been interpreted 
and/or predicted within the framework of MCT. Sim- 
ulation results in this article and in Ref4i by varying, 
respectively, the size disparity and mixture composi- 
tion, provide a global picture which might motivate 
experiments in a huge variety of systems in order to 
test eventual analogies and crossovers between differ- 
ent relaxation scenarios as those presented here. 
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